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ABSTRACT 

The existence of supermassive black holes (SMBHs) as early as z ~ 7 is one of the great unsolved 
problems in cosmological structure formation. One leading theory argues that they are born during 
catastrophic baryon collapse in z <~ 15 protogalaxies in strong Lyman- Werner (LW) UV backgrounds. 
Atomic line cooling in such galaxies fragments baryons into massive clumps that are thought to 
directly collapse to 10 4 - 10 5 Mq black holes. We have now discovered that some of these fragments 
can instead become supermassive stars that eventually explode as pair-instability supcrnovae (PI SNc) 
with energies of ~ 10 55 erg, the most energetic explosions in the universe. We have calculated light 
curves and spectra for supermassive Pop III PI SNe with the Los Alamos RAGE and SPECTRUM 
codes. We find that they will be visible in near infrared (NIR) all-sky surveys by WFIRST and WISH 
out to z ~ 20, perhaps revealing the birthplaces of the first quasars. 

Subject headings: black hole physics - cosmology: early universe - theory - galaxies: formation - 
galaxies: high-redshift - stars: early-type - supernovae: general - radiative transfer 
- hydrodynamics - shocks 



1. INTRODUCTION 

One model for the origin of SMBHs, which have now 
been found at z ~ 7, o r less than a Gyr after the big bang 
(|Mortlock et alj|2011l) , is catastrophic baryon collapse in 
protogalaxie s in s trong LW UV backgroun ds at z ~ 15 
(iWise et al.ll200l [Regan fc Haehnel3l200l IShang et al.1 
In this scenario, the primitive galaxy is built up 
by mergers between halos and by accretion in the vicin- 
ity of nearby LW UV sources that completely suppress 
star formation in the halq s without evaporating them 
(jWolcott-Green et al.l l20lTl ). When the galaxy reaches 
~ 10 s Mq, its virial temperature crosses the threshold 
for atomic hydrogen line emission and its baryons begin 
to rapidly cool and collapse. Infall rates at the center 
of the galaxy can be enormous: 0.1 - 1 Mq yr" 1 , or 
1000 times those in which the first stars form at z ~ 
25. Numerical simulations show that the baryons can 
shed angular momentum via the "bars within bars" in- 
stability on multiple spatial scales and collapse into an 
isothermal atomically-cooled disk. New simulations show 
that such disks can either feed a single massive cen- 
tral object or fragment into several slightly sm aller ones 
(|Regan fc Haehnelti 120091: IWhalen et al.ll2012e[) . In Fig- 
ure [T] we show the formation and fragmentation of such 
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a disk at the center of a 10 8 Mq halo in a strong LW 
background at z ~ 15 in an adaptive mesh refinement 
(AMR) simulation done with Enzc0. 

The evolution of the fragments depends on their 
masses at birth and subsequent accretion histories (see 
lOhkubo et al.l 120091 for studies of Pop III stellar evo- 
lution under ongoing accretion at much lower rates). 
One possibility is that a black hole forms at the cen- 
ter of the fragment early in the collapse that ra- 
diatively supports its surrounding envelope, a " qua- 
sistar" (iBegelman et all [20061 [20081: IBegelman! [20101: 
iVolonteri fc Begelmanl^OiaT The fragment may instead 
gradually collapse while being partially supported by un- 
stable nuclear burning, never reaching the main sequence 
or becoming a star. In this scenario, transient nuclear 
burning fails to offset gravity, and the fragment c ollapses 
to a 10 4 - 10 5 Mq black hole ([Fuller et al.lll986[) . If the 
clump has sufficient angular momentum, collapse pro- 
ceeds via the formation of a black hole accretion disk sys- 
tem. If the fragment does form a star, recent calculations 
have shown that its UV radiation will probably not be 
able to halt accretion and it may evolve un der heavy in- 
fall over its entire life ([Johnson et al.ll2012b[) . In contrast, 
lower-mass Pop II I stars usually disp e rse the baryons 
from their halos (IWhalen et all 120041: iKitavama et al ' 



2004 lAlvarez et alJl200fit lAbel et a. .1 120071 IWise fc Abell 
2008T ). We note that other mechanisms for massive frag- 



mentation and SMBH seed formation have been pro- 
posed, such as the s uppression of gas cooling by primor- 
dial ma gnetic fields (|Sethi et aLll2010f) and cold accretion 
shocks ([Inavoshi fc Omukaill2012l ). 

We have found that for a narrow range of mass 
around 55,000 M Q , atomically-cooled fragments can set- 
tle into stable nuclear burning and become supermas- 
sive stars with lifetimes of ~ 2 Myr. These stars 
die as extremely energetic PI SNc, with energies of 

9 http:/ /code. google. com/p/enzo/ 
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Fig. 1. — Baryon collapse at the center of a z ~ 15 protogalaxy in a LW UV background in the Enzo AMR code l lWhalen et al. 2012c). 
Left: the formation of an atomically-coolcd accretion disk at the center of the nascent galaxy at 10,000 yr. Right: the breakup of the disk 
into several supermassivc fragments shortly thereafter. The scale is ~ 2000 AU. 

- 10 55 erg or 100 times those of 65 - 260 M fl Pop 
III PI SNe (IRakavv fe Shavivl Il967t IBarkat et dj |l963 
IHeeer fc Wooslevl I2002t IGal-Yam et al.1 I2009D . Such 
events would be the most energetic explosions in the uni- 
verse, and their detection could reveal the birthplaces of 
SMBH created by direct collapse, since LW protogalaxies 
are the only environments known to form such massive 
fragments at z ~ 15. Could suc h SNe be discovered by 
existing or future observatories? iWhalen et al.l (|2012ctf ah 
recently found that 140 - 260 M@ Pop III PI SNe will be 
detected in the NIR out to z > 30 by the Jam es Webb 
Space Telescope ((JWST) IGardner et all [2006]) and to 
z ~ 15 - 20 in all-sky surveys by the Wide-Field Infrared 
Survey Telescope ( WFIRST) and the Wide-Field Imag- 
ing Surveyor for High Redshift ( WISH). However, super- 
massivc explosions might occur in extremely dense ac- 
cretion envelopes that quench their luminosities at early 
times when the greatest fraction of their flux is redshifted 
into the NIR. It is not clear if SNe in such dense shrouds 
could be detected. 

We present numerical simulations of light curves, spec- 
tra and NIR signals of supermassivc Pop III PI SNe at 7 
< z < 30 done with the Los Alamos RAGE and SPEC- 
TRUM codes. We consider only the observational sig- 
natures of these events, and defer detailed discussion of 
progenitor evolution, explosive nucleosynthesis and mul- 
tidimensional mixing to a forthcoming paper. In § 2 we 
review our numerical methods for evolving the star, its 
explosion, and the propagation of the blast through the 
star and its surrounding envelope. In § 3 we examine 
blast profiles, light curves and spectra for the SN in the 
source frame, and in § 4 we show NIR light curves in the 
observer frame and determine detection thresholds as a 
function of redshift for these explosions. In § 5 we con- 
clude by discussing complementary detection strategics 
for the direct collapse of supermassivc disk fragments and 
the formation of SMBH seeds. 

2. NUMERICAL METHOD 

We calculate light curves and spectra in three stages. 
First, we evolve the 55,500 M zero metallicity star from 



the beginning of the main sequence through core collapse, 
nuclear burning, and the arrival of the shock at the edge 
of the star in the Kepler code. In a parallel calculation 
we map a Kepler profile of the star onto a 2D AMR grid 
in the CASTRO code at the onset of collapse and evolve 
it through core bounce, explosive burning and again up 
to just before shock breakout. Second, we spherically av- 
erage mass fractions from the final CASTRO profile onto 
a ID AMR grid in RAGE together with final density, ve- 
locity and energy profiles from the Kepler calculation. 
This is done to model how mixing in the interior of the 
star prior to shock breakout alters explosion spectra at 
later times. We evolve the shock through the surface of 
the star and into the surrounding medium with RAGE 
until it dims below observability. Finally, we post pro- 
cess our RAGE profiles with the SPECTRUM code to 
calculate light curves and spectra. 

2.1. Kepler 

We calculate the internal structure of the star by evolv- 
ing it from the zero age main sequence to the begin- 
ning of collapse in the ID stellar evolutio n code Kepler 
(|Weaver et al.llT97l IWooslev et all[200l . The SN be- 
gins when the core of the star encounters the pair insta- 
bility during central He burning, begi ns to contract, and 
trigg e rs explosive O and Si burning dHeger fc Wooslevl 
200a iScannapieco et al.l 120051; iHeger fc Wooslevl 12010: 



Joggerst fc Whaled 120111: iKasen et alJ l20lH iPan et al l 
2012allbl K The non-rotating star, which is resolved with 



1148 mass zones, lives for 1.69 Myr and then dies as com- 
pact blue giant with a radius of 1.33 x 10 13 c m, sim ilar to 
those of the z-series stars in IWhalen et ail (I2012bn even 
though they are ~ 200 times more massive. The mass 
of the He core at the time of the explosion is 2.67 x 10 4 
Mq. This treatment is approximate because the super- 
massive star would likely accrete gas throughout its life, 
and such flows will affect its evolution. If accretion onto 
the sta r is governed by radiative feedback (pohnson et al.l 
l2012bf) . much of its evolution will proceed at a constant 
mass and our approach will be valid. The evolution of 
massive primordial clumps and supermassivc stars under 



3 



ongoing accretion will be the focus of a future study. 

2.2. CASTRO 

At the beginning of core collapse we map our Ke- 
pler profiles onto a 2 D axisymmetric grid in CASTRO 
(|Ahngren et alJfeOlOt) and then evolve the SN through 
core bounce and explosive burning, halting the advance 
of the shock when it reaches the edge of the star in mass 
coordinate. CASTRO (Compressible ASTROphysics) is 
a multidimensional Eulerian AMR code with a unsplit 
Godunov hydrodynamics solver. Energy production is 
calculated with a 19-isotope network up to the point 
of oxygen depletion in the core and with a 128-isotope 
quasi-equilibrium network thereafter. We evolve mass 
fractions for the 15 even- numbered elements predomi- 
nantly synthesized by PI SNc. Radiation transport is 
not required in these simulations because the mean free 
paths of photons prior to breakout are so short that they 
are simply advected through the star by the fluid flow. 
We include the contribution of photons to pressure in the 
equation of state. Our models include energy deposition 
due to radioactiv e decay of 56 Ni in the e .jecta as described 
by equation 4 in Uoggerst et al.l ()2010f ) although, as we 
discuss below, this explosion produces very little 56 Ni, 
unlike 140 - 260 M Pop III PI SNe. 

Mapping an explosion profile from a ID Lagrangian co- 
ordinate mesh in mass to a 2D mesh in space can lead to 
violation of conservation of mass and energy. Linear in- 
terpolations in radius can also fail to resolve key features 
of the original profile, such as the structure of the core of 
the star and its temperature profile. Failure to properly 
map temperature features can be especially problematic 
because nuclear burn rates are highly sensitive to them 
during the explosion. To avoid these difficulties, we port 
Kepler profiles to C ASTRO with the new conservative 
mapping scheme of IChen et all ([20TTf ). This approach 
conserves mass and energy while reproducing all the fea- 
tures of the original profile over a broad dynamical range 
in space. The CASTRO root grid is 256 2 with a reso- 
lution of 2.0 x 10 10 cm and up to two levels of AMR 
refinement (a factor of four increase in resolution). 

The star explodes with an energy of 7.74 xlO 54 erg. 
Explosive burning begins in the O and Si layers and is 
finished by ~ 100 s. The SN creates only trace amounts 
of 56 Ni, ~ 2.25 x 10~ 8 M Q , unlike 140 - 260 M Q Pop III 
PI SNe that synthesize up to 50 M of 56 Ni. The core of 
the SMS does not burn all the way to 56 Ni during the ex- 
plosion as in less massive PI SNc. 56 Ni is only formed at 
the edge of the He layer by the shock as it is launched out 
into the star. As we show in Figure [2j the interior of the 
star is heavily mixed by the time the shock reaches the 
surface, un like low-mass Pop III PI S Ne that exhibit lit- 
tle mixing (jJoggerst fc W halcn 20TlT) . The mixing is due 
to fluid instabilities that are seeded during core collapse 
and then amplified by explosive burning rather than the 
later formation of a reverse shock and the subsequent ap- 
pearance of Raylcigh- Taylor instabilities, like those that 
mix the interiors of 1 5 - 40 Mq Pop III core-collapse SNc 
(jJoggerst et al.|[2010f ). Pre-breakout mixing is important 
to explosion spectra because it determines the order in 
which emission and absorption lines appear in the spectra 
over time. Because elemental mass fractions are realis- 
tically distributed in radius and angle in CASTRO just 



before shock breakout when we spherically average them 
and map them into RAGE, we can capture the approx- 
imate effect of mixing on spectral lines in ID radiation 
hydrodynamical calculations. We stop the CASTRO run 
when the shock is no less than 100 photon mean free 
paths A p from the edge of the star: 



where Kxh is the opacity due to Thomson scattering from 
electrons (0.288 gm _1 cm 2 for primordial gas), and p is 
the density just beyond the shock inside the star. 

2.3. RAGE 

We follow the shock through the surface of the star 
and its envelope with the Los Alamos National Lab- 
oratory (LANL) rad iation hydrodynamics code RAGE 
(|Gittings et al.ll2008t) . RAGE (Radiation Adaptive Grid 
Eulerian) is a multidimensional AMR code that cou- 
ples second order conservative Godunov hydrodynam- 
ics to grey or multigroup flux-limited diffusion (FLD) 
to model strongly radiating flows. RAG E utilizes the 
LAN L OPLIB atomic opacity databas<£3 (|Magee et all 
[19951) and can evolve multimaterial flows with a variety of 
eq uations of s t ate (E OS). We employ the same physics as 
in iFrev et alJ (|2012f l: multispecies advection, grey FLD 
radiation transport with 2-tcmpcrature (2T) physics and 
LTE opacities, energy deposition from the radioactive 
decay of 56 Ni, and no self-gravity. Evolving radiation 
and matter temperatures separately in RAGE more ac- 
curately models shock breakout and its aftermath, when 
they can be out of equilibrium. We advect mass fractions 
for 15 elements, the even numbered elements predomi- 
nantly synthesized by PI SNe. 

We spherically average densities, velocities, specific in- 
ternal energies (erg gm _1 ), and species mass fractions 
from CASTRO onto a 50,000 zone ID spherical AMR 
mesh in RAGE. Since radiation energy densities are not 
explicitly evolved in Kepler, we initialize them in RAGE 
as 

E rad = aT\ (2) 

where a = 7.564 x 10~ 15 erg cm~ 3 K~ 4 is the radiation 
constant and T is the gas temperature. Also, because the 
gas energy in Kepler includes contributions from the ion- 
ization states of atoms, we construct the specific internal 
energy from T with 

Egas — CyT, (3) 

where Cy = 7.919 x 10 7 erg K _1 is the specific heat of 
the gas. 

At the beginning of the simulation we resolve the re- 
gion from the center of the grid to the edge of the shock 
in the velocity profile with 5000 zones on the root grid. 
We allow up to five levels of refinement during the ini- 
tial mapping of the profile onto the grid and throughout 
the simulation. Our grid resolves all important density 
and velocity features in the setup with a minimum of 
10 zones and up to 32 times more zones if needed. It 

10 http:/ /aphysics2/www. t4.lanl.gov/cgi-bin/opacity/tops. pi 
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Fig. 2. — Mixing in the Ni shell (left), O shell (center) and Si shell (right) just before shock breakout in CASTRO. The images are mass 
fractions. 
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Fig. 3. — Initial RAGE explosion profiles. Right: densities. Left: velocities 
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also ensures that the photosphere of the shock is always 
resolved; failure to do so can lead to underestimates of lu- 
minosity during post processing. The radius of the shock 
in our setup is 1.3 x 10 13 cm and our initial root grid 
has a resolution of 2.6 x 10 9 cm with an outer boundary 
at 1.3 x 10 14 cm. 

We set reflecting and outflow boundary conditions on 
the fluid and radiation variables at the inner and outer 
boundaries of the grid, respectively. At the beginning of 



the simulation, Courant times are short due to high tem- 
peratures, large velocities and small cell sizes. To mini- 
mize execution times and to accommodate the expansion 
of the SN, we periodically regrid the profiles onto a larger 
mesh as the explosion grows. Each time we regrid the SN 
we allocate 5000 zones out to either the edge of the shock 
(prc-brcakout) or the edge of the radiation front (post- 
breakout). In the latter case we take the radius at which 
the radiation temperature falls to the wind temperature 
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(0.01 eV) to be the edge of the front. The inner bound- 
ary is always at the origin and the outer boundary of the 
final, largest mesh in our simulations is 1.0 x 10 18 cm. 

2.4. Circumstellar Envelope 

We consider explosions in two kinds of envelope: low- 
mass outflows (SMS1) and massive inflows like those that 
grew the star to such large masses in such short times 
(SMS2). For diffuse outflows we adopt the usual power- 
law density profile for a wind at a constant velocity: 

Pw (r) = — — . 4 

Here, m is the mass loss rate associated with the wind 
and vw is the wind speed. The mass loss rate is cal- 
culated from the total mass loss M tot and the main se- 
quence lifetime of the star £msl: 



Since Pop III stars arc not thought to lose much mass 
over their lives because th ere are no metals in their at- 
mospheres to drive winds dKudrit zki 2000; Baraffe et al.l 
[20011 IVink et al.l [200ll IKrticka fc Kubatl 12006ft . we set 
M tot = 0.1 M G and v w = 1000 km s" 1 . 

We treat the massive infall envelope as a wind in re- 
verse, with m = 0.01 Mq yr^ 1 and an infall velocity 
v w = 5 km s _1 , in keeping with numerical simulations 
of baryon collapse in protogalaxies in strong LW back- 
grounds. This envelope assumes that accretion is spher- 
ical when in reality it may occur in a disk, so it should 
be considered to be the most massive one through which 
the SN shock might propagate. In both cases we take H 
and He mass fractions to be 76% and 24%, respectively, 
and for simplicity we take the wind to be cold (T = 0.01 
eV) and neutral in all our models. 

This latter assumption holds for dense envelopes, 
where semi-analytical studies have shown that ionizing 
UV photons cannot propagate more tha n a few dozen 
stellar radii from the star over its lifetime (jJohnson et al.l 
l2012b|) . On the other hand, the compact blue progenitor, 
with a total luminosity of 3.5 x 10 42 erg s _1 and T e g = 
6.85 x 10 4 K, certainly ionizes the diffuse wind, so the 
light curves and spectra we calculate for that case should 
be taken as lower limits. We show initial RAGE density 
and velocity profiles for the shock, the star, and its sur- 
rounding envelope in Figure |3] The surface of the star 
is evident as the sharp drop in density at ~ 1.4 x 10 13 
cm. The fact that the infall density is more than seven 
orders of magnitude greater at the surface of the star 
than the outflow density has important consequences for 
shock and radiation breakout, as we discuss below. In 
both cases we evolve the SN out to 3 yr. 

2.5. SPECTRUM 

We calculate spectra for the explosions with the LANL 
SPECTRUM code. SPECTRUM directly sums the lu- 
minosity of every fluid element in a SN profile to calcu- 
late the total flux escaping the ejecta along the line of 
sight for 14900 wavelengths. The procedure is described 



in detail in iFrev et all (|2012f ) and accounts for Dopplcr 
shifts and time dilation due to the relativistic expansion 
of the ejecta. SPECTRUM also calculates the intensities 
of emission lines and the attenuation of flux along the 
line of sight with OPLIB opacities, so it captures limb 
darkening and absorption lines imprinted on the flux by 
intervening mate rial in the SN ejec ta and envelope. 

As explained in lFrev et al.l (|2012f) . densities, velocities, 
radiation temperatures and mass fractions from the finest 
levels of refinement in the RAGE AMR hierarchy are ex- 
tracted and ordered by radius into separate files, with 
one variable per file. These profiles usually contain more 
than 50,000 radial zones, so limits on machine memory 
and time prevent us from using all zones to calculate a 
spectrum. We therefore map only a subset of the points 
onto the new grid. We first sample the radiation energy 
density profile inward from the outer boundary to find 
the position of the radiation front, which we define to be 
where aT 4 rises above 1.0 erg/cm 3 . This energy density 
is intermediate to that of the cold wind and the first pho- 
tons to reach a given radiusin the envelope. The radius 
of the r = 20 surface is then found by integrating the op- 
tical depth due to Thomson scattering inward from the 
outer boundary, where n^h = 0.288 gm" 1 cm 2 for pri- 
mordial H and He. This gives the greatest depth from 
which photons can escape from the ejecta because kt/i is 
the minimum total opacity. 

To compute a spectrum, we interpolate the densities, 
temperatures, velocities and mass fractions we extract 
from RAGE onto a 2D grid in r and /x = cos 9 in SPEC- 
TRUM whose inner and outer boundaries are and 10 18 
cm, respectively. The region from the center of the grid 
to the r = 20 surface is partitioned into one hundred 
uniform zones in log radius. We allocate 5000 uniform 
zones in radius between the t — 20 surface and the edge 
of the radiation front. The wind between the front and 
the outer boundary of the grid is divided into 100 uni- 
form zones in log radius, for a total of 5200 radial bins. 
The fluid variables in each of these new radial bins is 
mass-averaged to ensure that SPECTRUM captures very 
sharp features in the original RAGE profile. The grid is 
discretized into 160 uniform zones in \i = cos 9 from -1 to 
1. Our choice of mesh yields good convergence in spec- 
trum tests, fully resolving regions of the flow from which 
photons can escape the ejecta and only lightly sampling 
those from which they cannot. 

3. BLAST PROFILES, LIGHT CURVES AND SPECTRA 

We show velocity and gas temperature profiles at shock 
breakout for the SMS1 and SMS2 explosions in Figures 
|4] and [5l Before breakout, the SN cannot be seen by 
an external observer because photons from the shock are 
scattered by e~ in the upper layers of the star. When the 
shock reaches the surface of the star it abruptly acceler- 
ates, as shown in the velocity profiles of Figures 0] and [5j 
The shock also releases a brief, intense pulse of photons 
into the envelope. This transient, which is mostly x-rays 
and hard UV, blows off the outer layers of the star, which 
detach from and accelerate ahead of the shock as we show 
at 8190 and 8540 s in the SMS1 velocities and at 8720 
and 1.01c04 s in the SMS2 velocities. This effect is more 
pronounced in SMS1 because it is easier for the radia- 
tion front to drive a precursor into the diffuse wind than 
the dense infall. The advancing radiation front is visi- 
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Fig. 4. — Shock breakout into a diffuse envelope (SMS1). Right: velocities; from left to right: 7980 s, 8050 s, 8190 s and 8540 s. Left: 
temperatures at the same times. 
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Fig. 5.— Shock breakout into dense infall (SMS2). Right: velocities; from left to right: 7980 s, 8100 s, 8720 s and 1.01e4 s. Left: 
temperatures at the same times. 



ble as the flat plateau in gas temperature that extends 
from the outer edge of the shock into the surrounding 
medium. The temperature to which the radiation heats 
the gas falls as the shock expands, cools, and its spectrum 
softens (note that the temperature of the shock itself is 
much higher). 

At breakout there are marked differences in the profiles 
of the two explosions as the shock plows into the enve- 
lope. First, in SMS1 the shock accelerates as it exits the 
star and then continues to speed up as it descends the 
density gradient of the diffuse wind. In SMS2 the shock 
accelerates but then slows down as it crashes out into the 
much denser infall envelope. The radiation front also ad- 
vances more quickly into the diffuse wind than the infall. 
On the other hand, when the shock breaks out into the 
dense envelope it heats it to much higher temperatures. 
This hardens the spectrum of the shock and raises the 



temperature of the surrounding gas to higher tempera- 
tures than in SMS1, ~ 100 eV instead of ~ 50 eV. 

We show bolometric light curves for SMS1 and SMS2 
in Figure [5] The light curve for SMS1 is similar to those 
of lower-mass Pop III PI SNe except that the explosion 
is approximately ten times as luminous, as shown by the 
overlaid u250 and z250 light curves from IWhalen et al.l 
(|2012bf) . SMS1 exhibits the classic breakout transient, 
whose width is related to the light-crossing time of the 
star but is somewhat broader due to ra diation-matter 
couplin g effects as discussed in §4.1 of IWhalen et al.l 
(|2012bf) . Its light curve is similar in structure to that 
of z250, which is also the explosion of a compact blue 
giant (u250 is the PI SN of a 250 M Q red hypergiant). 
The double peak at early times common to all 150 - 250 
M Q z-series PI SNe and the 55,500 M Q PI SN is due to 
opacity effects. At breakout, the shock is extremely hot 
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Fig. 6. — Bolomctric luminosities for SMS1 and S MS2, together 
with lig ht curves for the u250 and z250 PI SNc from Whalen ct al. 
(20123). 
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Fig. 7. — Radiation breakout in SMS2: gas temperatures at 
3.69e6 s, 5.59e6 s, and 7.38e6 s. 

and the thermal transient ionizes the envelope. Later, as 
the fireball expands it cools, and the opacity of the shock 
and the envelope fall and allow more photons to stream 
out into space. The double peak is absent in the u-serics 
SNe because the shock is much cooler at breakout and 
never fully ionizes the envelope. At early times the lu- 
minosity of SMS1 comes from the conversion of kinetic 
energy into thermal energy by the shock. Later, its lu- 
minosity is due mostly to ejecta cooling, in contrast to 
low-mass PI SNe which are powered by 56 Ni decay. As 
in the z-series PI SNe, there is a resurgence in luminos- 
ity at ~ 10 7 s that is again due to opacity. At this time 
the t = 1 surface associated with the wavelength of peak 
emission in the spectrum has sunk to a hot layer deep in 
the ejecta, exposing it to the IGM and causing the SN 
to rebrighten. 



Radiation breakout in SMS2 occurs far after shock 
breakout, at ~ 5 x 10 6 s as we show in Figure [6] Radia- 
tion escapes the dense envelope much later because of its 
large optical depth, and when it happens it is gradual, as 
we show in Figure [7J Low-energy photons begin to leak 
out through the r = 1 surface for Thomson scattering at 
~ 2.5 x 10 16 cm by 3.7 x 10 6 s, and they are followed by 
more energetic photons by 5.6 x 10 6 s. At this point the 
shock is much cooler because of the large amount of PdV 
work it must perform on the dense shroud as it expands, 
but this results in an extremely luminous event in the 
NIR, as we discuss below. 

After radiation breakout the shock appears to flicker 
until ~2 x 10 7 s. This is due to radiative cooling and 
the cyclic formation and dissipation of a reverse shock 
in the ejecta. As the shock plows up the envelope a re- 
verse shock breaks free from the forward shock and is 
driven into the interior of the ejecta in the frame of the 
flow. As the reverse shock detaches and recedes from 
the forward shock, it loses pressure support to radia- 
tive cooling by emission lines in the shocked gas and 
retreats back toward the forward shock. As the for- 
ward shock continues to sweep up the envelope a reverse 
shock again forms and backstops from the forward shock. 
The cyclic heating and cooling of postshock gas asso- 
ciated with the oscillation of the reverse shock causes 
the fluctuations in luminosity from 5 x 10 6 - 2 x 10 7 s. 
The period of oscillation is determined by cooling rates 
in the gas (IChevalier fc Imamural I l982b llmamura et al.l 
ll984HAnninos et al.ll997|) and is independent of the mass 
swept up by the ejecta. Such ripples are also found in 
Lyman alpha emission by primordial SN remnants as 
they sweep up neutral ga s in cosmological hal os on larger 
scales (note Figure 11 in lWhalen et al.H2008f) . The light 
curves of both SMS1 and SMS2 are easily distinguished 
from those of less massive Pop III PI SNe. 

We show velocity and density profiles for both explo- 
sions at intermediate to late times in Figures [5] and |H1 
Multiple reverse shocks reverberate through the interior 
of SMS1 at earlier times and mostly dissipate by 3 yr, 
although some structures are still visible in the veloc- 
ity. These shocks are due to the internal dynamics of the 
ejecta rather than being driven inward by swept up gas 
(indeed, the remnant does not accumulate its own mass 
in ambient gas until it has grown to 6 pc). The density 
exhibits more self-similarity than the velocity, although 
its structure is always evolving in the vicinity of the shock 
from which most of the luminosity originates. The for- 
mation of a strong reverse shock due to plowcd-up gas 
can be seen in the SMS2 velocity profiles from 10 5 s to 
10 6 s, and it leads to multiple shocks at 10 7 s. This net- 
work of shocks is responsible for the oscillation in SMS2 
bolometric luminosities from 10 6 - 2 x 10 7 s in Figure HI 
The SMS2 density profiles are less self-similar over time, 
and there is significantly more evolution at the shock, 
as would be expected for ejecta plowing up a massive 
envelope. 

4. NIR LIGHT CURVES 

We calculate NIR light curves from our spectr a with 
the photometry code developed bv lSu et all (|2011l ). Each 
spectrum is rcdshiftcd prior to removing the flux that 
is absorbed by intervening neutra l hydrogen alo ng the 
line of sight using the method of iMadaui (|l995l ). The 
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Fig. 8. — Evolution of SMS1 at intermediate and late times. Right: 
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Fig. 9. — Evolution of SMS2 at intermediate and late times. Right: 
densities at the same times. 

spectrum is then dimmed by the required cosmological 
factors for a specified redshift. The least sampled data 
is linearly interpolated between the input spectrum and 
filter curve to model the light curve in a given filter. 

4.1. SMS1 

NIR luminosities are plotted for SMSf at z = 7, 10, 15, 
20 and 30 in the left panel of Figure HOI The SN will be 
visible to JWST at all epochs for over 1000 days and can 
be detected by WFIRST out to z > 20 for over 200 days. 
We include JWST for completeness, although given the 
much larger field of view of WFIRST and WISH and 
the low rates expected for this type of SN it is far more 
likely that such events will be detected in all-sky NIR 
surveys. SMS1 is extremely luminous in the NIR, with 
peak magnitudes ranging from 25 at z = 7 to 26.5 at z = 
20. These light curves exhibit far more variability than 
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velocities; from left to right: 10 s s, 10 6 s, 10 7 s and 10 s s. Left: 

their redshifted bolometric light curves might suggest, 
eliminating any possibility that these events would be 
mistaken for high-redshift protogalaxies. This variation 
is due to the expansion and cooling of the fireball in the 
source frame. 

The NIR profiles of SMS1 are easily disti nguished from 
those o f the u-series and z-series PI SNe in lWhalen et al.1 
(|2012b| ) at all redshifts (see Figures 9 and 10). The SMS1 
NIR light curves are similar in shape to those of z-series 
PI SNc, but z-scries NIR luminosities never reach magni- 
tudes brighter than 28 at z > 7. The SMS1 light curves 
evolve with redshift as expected: they broaden as z in- 
creases and the optimum filter wavelength increases with 
redshift. The NIR luminosities rise more quickly than 
they decline so these events are most easily detected in 
their early stages, but they nonetheless exhibit enough 
variability at any stage to be found in multiyear baseline 
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searches. 

4.2. SMS2 

We show NIR luminosities for SMS2 at z = 7, 10, 15, 
20 and 30 in the right panel of Figure [TU] They are quite 
different from those of SMS1. Consistent with radiation 
breakout from the shroud at ~ 58 days, no NIR signal 
is observed from these events until ~ 500 days at z > 
7. This explosion eventually becomes hundreds of times 
brighter in the NIR than SMS1, with peak AB magni- 
tudes from 21 at z = 7 to 23 at z = 15. It is visible 
to JWST for 200 days at all redshifts and for 200 days 
out to z > 20 by WFIRST and WISH. The much higher 
NIR luminosities are due both to the large radius of the 
fireball at radiation breakout and the relatively low tem- 
perature of the shock at this radius (~ 10 eV) because 
of the PdV work the fireball must do against the dense 
envelope. These lower temperatures drive the redshifted 
peak of the shock's spectrum into the NIR in the source 
frame. 

From 7 < z < 15, the steady rise in luminosity from 
500 - 700 days is due to the expansion of the surface area 
of the fireball, and the relative magnitudes of the three 
light curves are properly ordered in redshift. The shal- 
low ripples in luminosity over these redshifts have much 
shorter periods than those in the bolometric luminosity 
in Figure [6] and are likely due to opacity fluctuations in 
the shock. The sudden drops in luminosity at 675 days 
and 700 days are due to sudden jumps in opacity that 
are associated with the cooling of the fireball. At z > 15, 
the SN is actually brighter in the JWST Mid-Infrared In- 
strument (MIRI) bands after 700 days than in the NIR, 
as we show in Figure 1111 

In sum, SMS explosions in both diffuse winds and 
dense envelopes will be visible in JWST NIR deep fields 
out to z > 30 and in all-sky NIR surveys by WFIRST 
and WISH out to z ~ 20, which is earlier than their 
likely era of formation. The fact that these NIR profiles 
change considerably with circumstellar envelope suggests 
that they will be powerful probes of the environments 
of such explosions. It is worth noting that even fully 
shrouded explosions will be visible at the earliest epochs. 
The envelopes we have chosen should bracket those in 
which SMS PI SNc will occur, so the NIR signals of ac- 
tual explosions may be intermediate to those of these two 
events. As noted above, both SMS1 and SMS2 are eas- 
ily distinguished fro m low-mass Pop III P I SNe as well 
as co re-collapse SNe (|Whalen et al.ll2012dD an d Type Iln 
SNe (IMoriva et alJl2010b IWhalen et al.H2012ah . 

5. CONCLUSION 

The discovery of supcrmassivc Pop III PI SNe would 
confirm for the first time that massive fragments capa- 
ble of collapsing to 10 4 - 10 5 M SMBH seeds do form 
in primeval galaxies at high redshift. Although the rate 
of such events remains unknown, it might be thought 
that they are very rare because supermassive fragments 
must fall into a relatively narrow mass range to actu- 
ally become stars and because few protogalaxies form 
in LW backgrounds capable of fully sterilizing them of 
H 2 . However, recent developments suggest that these 
processes were more frequent than previously thought. 

First, new simulations indicate that the assembly of 
protogalaxies in strong LW backgrounds may have been 



relatively common, yielding higher rates of SMBH seed 
production than might naively be inferred from the num- 
ber d ensity of z ~ 7 quasars, ~ 1 Gpc -3 ([Agarwal et al.l 
2012). The sustained exponential growth required to 
reach such masses depe nded on the topology of cold 
flows over cosmic time (|Di Matteo et al.l 120121 ) . so the 
scarcity of such flows may have governed the density of 
z ~ 7 quasars, not the rate of seed formation. Second, 
rotation could broaden the mass range over which su- 
permassive stars encounter the pair instability by en- 
hancing mixing and more rapidly build ing up massive 
He cores (jChatzopoulos &: W heeler 2012). Greater mass 
ranges imply larger event rates. Indeed, rotation, whose 
effects were not considered in the first studies of super- 
massive stellar evolution or in our simulations, may also 
increase the mass range over which supermassive frag- 
ments can actually become stars by supporting more of 
them against collapse and enabling stable nuclear burn- 
ing. 

A reasonable upper limit to SMS PI SN event 
rates are those of 1 40 - 26 M Q Pop III PI SNe, 
which IWhalen et al.l (|2012bl ) and others find to be 
~ 10~ 2 yr _1 deg -2 at z > 10, which implies 
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be well below this limit, precluding their detection 
by JWST, they will clearly be bright enough to ap- 
pear in wide field campaigns. They might also be de- 
tected in present NIR all-sky surveys by the Subaru 
Hyper Suprime-Cam at z > 3 but furthe r calculations 
will be necessary to con firm this (e.g., iTanaka et all 
120121 IMoriva et all [2012h . Synchrotron emission from 
Pop III SNe at z > 10 can be detected at 21 cm 
by existing observatories such as eVLA and eMerlin 
and fut ure ones such as the Sq uare Kilometer Array 
(SKA) (jMeiksin fc Whalenl \20m . SMS PI SN explo- 
sions in dense envelopes may likewise generate enough 
synchrotron emission to be discovered in radio surveys. 

Given that most supermassive stars will still directly 
collapse to BH, could there be other ways of detect- 
ing SMBH seed formation in protogalaxies? Past stud- 
ies have shown that collapsing supermassive stars be- 
come extremely luminous in thermal neutrino emission 
as the central BH forms, w i th en ergies of ~ 10 MeV 
([Fuller et al.l [l98l iShi et all Il998l) . The prospects for 
detecting such neutrinos depends on the initial mass and 
entropy of the core. Although the total energy emitted 
from these massive stars increases linearly with mass, the 
entropy of their cores also increases with mass. Higher 
entropies lead to larger proto-blac k holes with lower peak 
densities and lower temperatures. iFrver fc Hegerl ([201 ID 
found that the neutrino luminosity does not increase 
much with mass for stars above 10,000 Mq. The mean 
electron neutrino energy for stars above 10,000 Mq is 
below 6 MeV and the \x T energy is not much higher. The 
collapse of such cores would be difficult to observe with 
neutrino detectors. However, if the core entropy is closer 
to that of a 1,000 Mq star, the luminosity peaks more 
dramatically. Although even for these cores the mean 
electron neutrino energy is ~ 7 — 8 MeV, the fi T energy 
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Fig. 10. — JWST NIRCam light curves for supcrmassivc Pop III SNc in diffuse winds (SMS1, left panel) and dense envelopes (SMS2, 
right panel). The horizontal dashed lines at AB magnitudes 32 and 27 are the photometry limits for JWST and WFIRST, respectively. 
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ISurman et all 120061 |20"o7I). The nucleosynthetic signa- 
ture of this process, which could be imprinted on stars 
that later form in the protogalaxy, depends on the tem- 
perature of the disk and hence the radius of the BH. 
It can therefore provide a diagnostic of the mass of the 
SMBH seed at birth, since massive BH with large event 
horizons burn at lower temperatures and yield chemi- 
cal abundances that are distinct from those of smaller 
disks, which can burn all the way to Ni. Ancient, dim 
metal-poor stars bearing the ashes of this process could 
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Fig. 11.— JWST MIRI NIR light curves for SMS2. 



lies in the 20-30 MeV range and would be more easily de- 
tected after cosmological rcdshifting. More detailed cal- 
culations are needed to be certain, but these cores would 
likely contribute to the neutrino background in detectors 
such as IceCube. If the density profile of the collapsing 
star also imposes a unique spectrum on the neutrino flux, 
it would facilitate its extraction from noise. 

As noted in the Introduction, collapse may also lead 
to the formation of a black hole accretion disk system, 
with nuclear burning near the event horizon whose prod- 
ucts could be expelled out into the halo by a jet (e.g., 



supcrmassive star could also emit gravity waves (GWs) 
that migh t be detected in existing or future G W facil- 
ities (e.g.. iFrver eTail [20021: IFrver fc New|[20ll . These 
multi-messenger events, together with the most energetic 
supcrnovae in the universe, may soon herald the births 
of the first quasars. 
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